A Phase Glass is a Bose Metal: New Conducting State in Two Dimensions 
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In the quantum rotor model with random exchange interactions having a non-zero mean, three 
phases, a 1) phase (Bose) glass, 2) superfluid, and 3) Mott insulator, meet at a bi-critical point. We 
demonstrate that proximity to the bi-critical point and the coupling between the energy landscape and 
the dissipative degrees of freedom of the phase glass lead to a metallic state at T = 0. Consequently, 
the phase glass is unique in that it represents a concrete example of a metallic state that is mediated 
by disorder, even in 2D. We propose that the experimentally observed metallic phase which intervenes 
between the insulator and the superconductor in a wide range of thin films is in actuality a phase glass. 
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There is now a preponderance of experimental evidence 
that the disorder or magnetic field-induced destruc- 
tion of superconductivity in a wide range of thin metal 
alloy films leads first to a metallic state with a non-zero 
conductivity as T — + 0. At sufficiently large values of 
the disorder or magnetic field, a transition to a true in- 
sulating state obtains. Within the standard bosonic de- 
scription |^|,|| of the insulator-superconductor transition 
(1ST), the onset of an intervening metallic state is prob- 
lematic because only two options are thought to exist 
for bosons: 1) localized in a Mott insulating state or 2) 
condensed in a superfluid. In the former, the conduc- 
tivity vanishes whereas the latter exhibits resistanceless 
transport. Further, including degrees of freedom which 
lie outside the bosonic or phase-only models, for exam- 
ple electronic excitations, is of no help as electrons are 
localized in 2D. Indeed, while the onset of the insula- 
tor in homogeneously disordered thin films is consistent 
with the emergence of electronic excitations, the in- 
tervening metallic and the subsequent superconducting 
states appear to be inherently bosonic in origin. 

Consequently, recent theoretical effort |]Io|-|l2| on the 
origin of the metallic state has focused strictly on bosonic 
models. Along these lines, we have shown [[ll] that 
the standard Mott insulating phase in a clean array 
of Josephson junctions has a non-zero conductivity as 
T — > 0. This result arises from the non-commutativity 
p3| of the frequency and temperature tending to zero 
limits of the conductivity in the vicinity of a quantum 
critical point, with u) = 0,T — > being the experimen- 
tally relevant limit for the dc conductivity. In the Mott 
insulator, quasiparticle excitations are gapped and obey 
a Boltzmann distribution. However, the collision time of 
such quasiparticles grows exponentially with the gap p| . 
As the conductivity is a product of the collision time and 
the quasiparticle density, the conductivity is necessarily 
finite in the limit lj < T . This type of Bose metal is 
fragile pd| ], however, and suppressed by dissipation and 
disorder. In fact, in the presence of disorder, the nature of 
the superfiuid-insulator transition changes dramatically. 
For example, several ||[l4| have argued that in the pres- 



ence of on-site disorder, destruction of the superfluid may 
(in the presence of incommensuration) obtain through an 
intervening phase with gapless excitations referred to as 
a Bose glass. In analogy with the Fermi glass, Fisher, et. 
al. |g| proposed that the Bose glass is an insulator with 
variable-range hopping conductivity. 

Nonetheless, we show explicitly that the glass phase 
which may interrupt the direct transition from a super- 
fluid to a Mott insulator in the generally disordered case 
is in fact a metal that has a well-defined T — > limit for 
the conductivity. As the thermal average of the super- 
conducting order parameter is non-zero but vanishes once 
averaged over disorder, we refer to the glass as a phase 
glass. We propose that the intervening metallic phase 
seen in the experiments is a phase glass. Our proof that 
such a phase possesses a non-zero conductivity as T — > 
constitutes the first demonstration of a stable metallic 
state in 2D in the presence of disorder. 

The starting point for our analysis is the charging 
model for an array of superconducting islands 
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with random Josephson couplings but fixed on-site 
energies, E@. The phase of each island is 6^. Note that 
additional on-site disorder of the form ivjd/dOj results in 
the equivalent particle-hole symmetric field theory pro- 
vided that the distribution of on-site energies has zero 
mean. The non-zero mean case is irrelevant here as this 
corresponds to a density-driven 1ST Jl5[ . Hence, our con- 
clusions apply to the general disordered case. To incor- 
porate ordered phases, we assume that the Josephson 
energies are random and characterized by a Gaussian dis- 
tribution, P{J i3 ) = l/V27rJ 2 exp -(Jij - J a ) 2 /2J 2 , with 
non-zero mean, Jo. The negative Josephson couplings 
included in this distribution are essential to the physics 
of a disordered superconductor |T^ ], particularly glassy 
ordering. We have studied the non-zero mean problem 
extensively and established explicitly the existence of 
a bi-critical point in which three phases meet, a Mott in- 
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sulator, phase glass, and superconductor. To distinguish 
between the phases, it is expedient to introduce the set of 
variables S, = (cos 0j, sin 6i) which allows us to recast the 
interaction term in the random Josephson Hamiltonian 
as a spin problem with random magnetic interactions, 
j) JijSi ■ Sj. Let (...) and [...] represent averages over 
the thermal degrees of freedom and over the disorder, re- 
spectively. In the superconductor not only (Sw) but also 
[(£«/)] acquire a non-zero value. In the phase (or spin) 
glass, however, (Si V ) ^ but [(£«,)] = 0, whereas in the 
paramagnet or Mott insulator, (Sw) = 0. 

The Landau theory for this problem is easily obtained 
|l0f using replicas to average over the disorder and the 
identity ln[Z] = Urri n _»o([i? n ] — l)/n to obtain the zero 
replica limit. The quartic and quadratic spin-spin inter- 
action terms that arise from the disorder average can be 
decoupled by introducing the auxiliary real fields, 



Q2t(k,wi,W2) = f3(2ir) d 5 d (k)5^ [£>M<W"2,o<J«6 
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Q«t(k,k',r,r') = (^(k,r)^(k',rO) 



(2) 



and 



*°(k,r) = (S*(k,r)), respectively. The super- 
scripts represent the replica indices. A non-zero value of 
^"(k, r) implies phase ordering in a charge 2e conden- 
sate. Hence, ^ couples to the charge degrees of freedom. 
For quantum spin glasses, it is the diagonal elements of 
the Q-matrix D(t - r') = lim, w0 jk (Q°°(k, k', r, r')) 
in the limit that |r — t'| — > oo that serve as the effective 
Edwards- Anderson spin-glass order parameter 16 lSy 
within Landau theory. The free energy per replica 

^[*,Q] = ^ SG (g)+ (fc 2 +^ + m 2 )i^(k, w „)i 2 

~ l d d X J d Tl dT 2 J2 ^(x,T 1 )^l(x,T 2 )Qf v (x > T 1 ,T 2 ) 

+U j drYsKM^T)] 1 (3) 

a, ii 

consists of a spin-glass part which is a third-order func- 
tional of the Q— matrices discussed previously ]l0| , |l6"| | , 
the terms that describe the charge 2e condensate and 
the term which couples the charge and glassy degrees 
of freedom. The parameters, k, t and U are the stan- 
dard coupling constants in a Landau theory and m 2 is 
the inverse correlation length. In the disordered phases, 
(\&u) = 0; hence, in the glassy phase, it is the fluctua- 
tions of the fy® field that survive. Our previous analysis 
|jl0f shows that the cross term we have retained here is 
the most dominant of the possible coupling terms near 
the bi-critical point. 

Our goal now is to calculate the charge transport in the 
glassy phase. Near the spin-glass/superconductor bound- 
ary m 2 should be regarded as the smallest parameter. 
Hence, it is the fluctuations of rather than those of 
Q ab that dominate. Consequently, we adopt the most 
general mean-field ansatz MJul for the Q— matrices 



The diagonal elements of the Q-matrices describe the ex- 
citation spectrum. In the glassy phase, the spectrum is 
ungapped and given by D(u>) = — \lo\/k. The linear de- 
pendence on \u>\ arises because the correlation function 
Quu(t) decays as t~ 2 |^6| , |l7| . This dependence results 
in a fundamental change in the dynamical critical expo- 
nent from z = 1 to z = 2 and the onset of overdamped 
dynamics (zero energy modes), thereby eliminating the 
5(lu) term from the conductivity. In addition, without 

loss of generality, we work in the replica symmetric case, 
q ab = q 

for all a and b as it was shown fl(],[l6| that replica 
symmetry breaking vanishes as T — > and our emphasis 
is the low-temperature limit. Finally, because our focus 
is charge transport and the electromagnetic gauge cou- 
ples only to the \E'° field, we retain only those terms in 
the free energy in which at least one of the ^ fields 
is present. Substituting the Q-matrix ansatz (Eq. (^)) 
into Eq. (^) and introducing a one-component complex 
field ip a — (^1,^2)1 we arrive at the following Gaussian 
theory, 

^"gauss= (k 2 +ojl + v\un\+m 2 )\tf; a (k,uj n )\ 2 

-0q Y *«n,0^ o Ck,w w )[^(k,w«)]* J (5) 

In the above action we introduced the effective dissipa- 
tion r] = l/(K 2 t) and rescaled q — > qnt. The associated 
Gaussian propagator is 

G^(k,u> n ) = G (k,u n )8 ab + /3G 2 (k,^„)g<L„,o (6) 

in the n — > limit Q with G (k,w n ) = (k 2 + uj 2 n + 
?7|aj| + m 2 ) -1 . The first term in Eq. (|J) is the stan- 
dard Gaussian propagator in the presence of Ohmic dis- 
sipation. The Ohmic dissipative term in the free-energy 
arises from the diagonal elements of the Q— matrices. 
However, it is the q— dependent term in the Gaussian free 
energy, the last term in Eq. (^|), that is new and changes 
fundamentally the form of the propagator. Because of 
the £ w „,o factor in the second term in the free energy, the 
propagator now contains a frequency-independent part, 
/3Gg(k, ui n = Q)q. In the free energy, this term couples 
different components of the replicas and hence cannot be 
regrouped with the mass term, m 2 . In fact, this term is 
a highly relevant perturbation in all dimensions. From 
simple tree-level scaling, ip — M d+2+2 )/ 2 -0', fc = k'/b,uj — 
ci//6 z , wc find that the term proportional to q in the 
free energy rescales as q' — qb 2+z . The dynamical ex- 
ponent z is determined by the fact the scaling dimen- 
sion of 77 should remain unchanged. This gives at the 
tree level z = 2 as proposed previously J(| for the su- 
perfmid/Bose glass transition. Hence, q' = qb 4 , implying 
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that the coupling to the energy landscape of the phase 
glass is strongly relevant and ultimately responsible for 
the metallic phase. 

To see how this comes about, we use the generalization 
p0| of the Kubo formula for the replicated action and 
write the conductivity to one-loop order per replica in 
the Gaussian approximation as 
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The conductivity contains two types of terms. All terms 
not proportional to q have been evaluated previously H] 
and vanish as T — > 0. The terms proportional to q 2 van- 
ish in the limit n — ► 0. The only terms remaining are 
proportional to q and yield after an appropriate integra- 
tion by parts 
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k 2 x G 2 Q (k,0) [G o (k,0)-Go(k,w n )] 



The momentum integrations are straightforward and 
yield 



a(uj = 0,T^ 0) = 



<1>1 



h 2m 4 
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a temperature-independent value for the conductivity as 
T — > 0. The dependence on q and r\ implies that dissi- 
pation alone is insufficient to generate a metallic state. 
What seems to be the case is that a bosonic excitation 
moving in a dissipative environment in which many false 
minima exist does not localize because it takes an expo- 
nentially long amount of time to find the ground state. 
This is the physical mechanism that defeats localization 
in a glassy phase. Further, the conductivity scales as 
1/m 4 and hence diverges as the superconducting phase 
is approached. This is precisely what is seen experimen- 
tally §%. 

That the conductivity plateaus in the phase glass 
regime does not appear to have been anticipated previ- 
ously. We now appeal to much more general arguments 
to prove that the singular dependence of the conductivity 
on m 2 as T — > survives even in the presence of the quar- 
tic interaction. At the tree level, a dynamical exponent 
of z = 2 renders the quartic interaction U marginally 
irrelevant. However, considering the last term in Eq. 
(|J) on equal footing with U in the one-loop renormal- 
ization group scheme, we reach the conclusion that the 
RG equations flow to strong coupling. The relevance of 
q at all dimensions manifests itself also by the increasing 
singularity of relevant contributions from higher order di- 
agrams in the perturbation series in U. We consider first 
the linear U correction. At this level, the self energy is 
given by a standard tad-pole diagram that arises from the 



couplings in the average {il) a ip ^2 c ip c il)* c ip c il)* c ) 1 yield- 
ing S = U / w / k G^(k,u; n ). This diagram @ is regu- 
larizable only once the term w 2 is retained in the prop- 
agator. The first term in Eq. (^) leads at T — to a 
standard mass renormalization and innocuous logarith- 
mic corrections, while the last term is more singular, giv- 
ing S' 1 ' = Uq/{Airm 2 ). A similar analysis can be un- 
dertaken to find the first-order correction to the T = 
conductivity. The relevant diagrams psfl can be readily 
generalized for a two-replica propagator given by Eq. (|^). 
The straightforward evaluation of the contribution that 
does not vanish in the T — > limit yields 



' (1) (- = o) = lr^, 

v ' Ah TTUl 



(9) 



suggesting that each subsequent order in the interac- 
tion leads to a more singular contribution to the self- 
energy. This points to a scaling function of the form 
a w (e 2 /h)(?]q/m 4 )& (g/m 2 ) where $(y) ~ y p for large 
y, which yields the critical behavior a ~ m~ x with 
x = 4 + 2p. The value of the exponent p cannot inferred 
at any finite order in perturbation theory. 
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FIG. 1. Diagrammatic representation of the conductivity. 
Each solid line denotes the fully renormalized propagator, G a t 
(see Eq.@) while the shaded rectangle is the vertex func- 
tion. The letters a, b, c represent the replica indices, and the 
internal momenta ki are not shown for simplicity. 



Nonetheless, we assume that all of the most singular 
diagrams can be resummed. A simple inspection of the 
perturbation series suggests that the fully renormalized 
propagator, 



G ab (k, u n ) = G(k, uj n )Sab + (3qg(k.)6 Uni o, 



(10) 



can be broken into replica diagonal and off-diagonal 
pieces. Likewise, we define the self-energy associated 
with this propagator to be S a b(k, ui n ) = E(k, w„)<5 a h + 
Pq9Qs.)S Uni o which contains formally all interaction 
terms. From the Dyson equation G a b — G^ + 
Gac^cdGdbj hi which the summation over the repeated 
indices is implied, we have that G(k,o; n ) = (G " 1 (k,w„) — 
f^k,^))- 1 and 5 (k) = [l + 0(k)]/(G o - 1 (k,O)-E(k,O)) 2 . 
The renormalization of the interaction U leads to the 
appearance of the corresponding vertex function p£(| 
r(ki, k2, k; fii, O2, w n ) which is connected to the self- 
energy by means of the standard Dyson equation. This 
vertex function enters the general expression for the con- 
ductivity represented diagrammatically in Fig. (nh. Each 
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solid line represents the renormalized propagator, G a &, 
while the shaded region denotes the vertex function, 
r(ki, k 2 , k = 0; Qi,D,2, w„). We are interested here only 
in the static T = conductivity. Once the first term 
in this diagrammatic expansion is integrated by parts, 
use of the standard Ward identity leads immediately to 
a cancellation of all diagrammatic contributions to a in 
which the external frequency vanishes. As a consequence, 
we obtain the leading contribution to the conductivity in 
the limit that u> = 0, T — > simply from a Taylor expan- 
sion around lj = 0. Using Eq. ( |l0| ) for the renormalized 
Green function, we obtain the exact expression, 
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+G(k 2 ,0) — (k l5 k 2 , |OD||o=o 



(11) 



for the temperature-independent part of the conductiv- 
ity. Here f (ki,k 2 , |w„|) = r(ki, k 2 , 0; -w n , u> n , w n ) + 
r(ki, k 2 , 0; 0, 0, uj n ), and it is taken into account that the 
frequency dependence of all functions enters through \u> n \ 
due to the full particle-hole symmetry. In deriving Eq. 
(pi]), we assumed that, 1) the infinite perturbation series 
in U is resumable in principle, and 2) that all propagators 
and the vertex function are analytic in \u> n \. The latter 
assumption seems reasonable, because the most singular 
contributions come from diagrams that do not contain 
frequencies at all. 

We have demonstrated here that the sluggish phase 
dynamics in a phase glass leads ultimately to a metal- 
lic state in d = 2 for bosonic excitations. The strong 
divergence of the resultant conductivity on m 2 is consis- 
tent with the experiments that have observed a distinct 
plateauing of the resistivity at low temperatures which 
increases in magnitude jjj as the distance from the true 
superconducting phase is increased. The metallic T = 
behavior obtains as a result of the coupling between the 
dissipative environment and the energy landscape of the 
phase glass. Further, since the dissipation inherent in a 
phase glass is independent of temperature, external dissi- 
pation arising from phonons is irrelevant as such coupling 
vanishes as T — > 0. Consequently, the metallic phase we 
have found here is robust to disorder and phonon scatter- 
ing and in fact constitutes the first explicit demonstra- 
tion of a metallic state in d = 2. Nonetheless, the theory 
presented here does not address the issue of whether or 
not the destruction of the superconducting phase occurs 
directly into a conventional insulator or a glassy phase. 
However, we suggest that only in the second scenario is 



the destruction of a 2D superconductor in the absence 
of a magnetic field consistent with the robustness of the 
metallic phase with respect to increasing disorder. We 
propose that aging and noise measurements as well as ex- 
periments sensitive to trapped flux should be performed 
in the intervening metallic regime to explore the glassy 
scenario suggested here. Clearly a promising extension 
of this work would be the fermionic case. 
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